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Abstract. We introduce an operational discord-type measure for quantifying nonclassical 
correlations in bipartite Gaussian states based on using Gaussian measurements. We refer to 
this measure as operational Gaussian discord (OGD). It is defined as the difference between the 
entropies of two conditional probability distributions associated to one subsystem, which are 
obtained by performing optimal local and joint Gaussian measurements. We demonstrate the 
operational significance of this measure in terms of a Gaussian quantum protocol for extracting 
maximal information about an encoded classical signal. As examples, we calculate OGD for 
several Gaussian states in the standard form. 
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1. Introduction 

Characterization and quantification of correlations in quantum systems are central in 
implementing quantum information processing tasks that cannot be done classically. Quantum 
discord was proposed as a measure of nonclassical correlations, which can capture 
correlations beyond quantum entanglement [HllHlSl. This measure of correlation was shown 
to be useful to characterize resources in a quantum computational model (DQCl) IH, quantum 
state merging nil a, remote state preparation [[71, encoding information onto a quantum 
state iSl, quantum phase estimation [[91, and quantum key distribution [fTOl . It was also shown 
that quantum discord is linked to entanglement generated by the activation protocol ifTTl or 
by a measurement [[T2l . 

In general, measures of quantum correlations can be defined as the difference between a 
quantum entropic measure and a classical entropic measure that is obtained from outcome 
probabilities of local measurements lfT3l . For a bipartite system in quantum state pab, 
quantum discord from subsystem i? to A is defined as 

D{B ^ A) = H^^^\A\B)-S{A\B), (1) 

where H^^''^\A\B) is the minimized classical conditional entropy, and the quantum 
conditional entropy is S{A\B) = S{A,B) — S{B), with S{A,B) and S{B) being the 

von Neumann entropies of the joint state and the marginal state pb = TI^a[pab]- The 

classical conditional entropy H^^^'^\A\B) is obtained from outcome probability distributions 
of measurements on A in the eigenbasis of the conditional states pA\b = Tl^B[pAB^b\/Pb 
obtained after performing a measurement on B described by POVM elements {!!(,}, where 
Pb = TrAB[pAB^b\ is the probability for outcome b [fTSl . Hence, we have 

= min (2) 

{nil ^ 

where 5'(y4|6) is the von Neumann entropy of pA\b and the remaining minimization is over the 
local measurements on B. In general, it is not clear how to perform this minimization for an 
arbitrary quantum state; however, it can be done for certain cases, including a large class of 
two-qubit states IfTTl . 

Quantum discord was generalized to quantify nonclassical correlations in continuous- 
variable systems, particularly Gaussian states. Interestingly, all Gaussian states except 
product states have nonzero quantum discord ifTSl [T^ . Gaussian quantum discord (GQD) 
was introduced as a measure of nonclassical correlations for Gaussian states in which the 
minimization in the classical conditional entropy Q is restricted to Gaussian measurements 
on B IfTTl [T^ . To get to the form Q, however, already requires nonGaussian measurements 
on the conditional states of A; these are generally measurements in a displaced squeezed 
number basis. Thus, the GQD cannot really be used as a figure of merit for Gaussian quantum 
protocols that only involve Gaussian states, Gaussian operations, and Gaussian measurements, 
such as a Gaussian version of the protocol in [|8][. 

By using the optimality of input Gaussian states for Gaussian channels lfT9l l20l . it was 
recently shown IfTTl that for a large class of Gaussian states, no nonGaussian measurements 
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on B can further minimize the value of quantum diseord, implying that Gaussian quantum 
diseord is equal to quantum diseord. It seems to be an open question whether this is true for 
all Gaussian states. 

In this paper, we introduee a new measure for quantifying nonelassieal correlations in 
bipartite Gaussian states, based solely on Gaussian measurements, which has qualitatively 
different behavior from GQD. This measure is defined as the differenee between the Gaussian 
version of the elassieal eonditional entropy H^^\A\B), whieh is given by minimizing over 
local Gaussian measurements on both subsystems, and the minimum conditional entropy that 
can be measured by a joint Gaussian measurement, H^qP\A\B). We refer to this measure 
as Operational Gaussian Diseord (OGD) because, firstly, it only depends on quantities that 
ean be measured via Gaussian operations and, seeondly, it has an operational signifieanee in 
terms of a quantum protoeol that is a Gaussian version of the protoeol in [l^. In this protoeol, 
a elassieal signal with a Gaussian probability distribution is encoded on one subsystem of a 
bipartite Gaussian state; using a local Gaussian or a joint Gaussian measurement, one tries to 
retrieve the signal from the noise assoeiated with the joint state and the measurement. The 
optimal measurement is the one that maximizes the elassieal mutual information between 
the measurement outeome and the input signal. We show that in the limit of large varianees 
for the signal probability distribution, the difference between the maximal elassieal mutual 
informations obtained by optimal joint and local Gaussian measurements is equal to the OGD 
of the bipartite Gaussian state. 

This paper is struetured as follows. In the following seetion, we review Gaussian 
quantum diseord. We introduee OGD in See. and ealeulate it for several Gaussian states 
in the standard form in See. We demonstrate the operational signifieanee of our measure 
in See. We eonelude the paper in the last section and pose an open problem. This paper is 
supplemented with one appendix. 

2. Gaussian quantum discord 

Quantum states that have Gaussian Wigner funetions are known as Gaussian states. We gather 
the phase-spaee quadratures for a two-mode system into a vector X = {xa,Pa, xb,Pb)- A 
Gaussian state pab can be fully characterized in terms of the mean quadratures {X) = X 
and the eovarianee matrix, whieh has elements [aABpj = ^{^iXj + XjXi) - (X,)(X,). 
As eorrelations in Gaussian states are independent of loeal displaeement operations, we ean, 
without loss of generality, assume that X = 0. Also, by applying local phase shifts and 
squeezing operations, any eovarianee matrix. 



( 3 ) 
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can be brought to the standard form in which A = diag(a,a), B = diag(6, 6), and 
C = diag(c, d) [|2l|23l, i.e., 


AB — 


/a 0 c o\ 
0 a 0 d 
c 0 6 0 
\0 d 0 bJ 


(4) 


Matrices A and B are the covariance matrices of the marginal states pA = TrB[pyis] and 
Pb = Tia[pab\, and the matrix C contains the quadrature correlations between the modes. A 
Gaussian state has zero discord if and only if C = 0 IfTSl . 

A closed-form expression for calculating the classical conditional entropy 
with restriction to Gaussian measurements on B was given in Ref. lITTl . A local Gaussian 
measurement is described by POVM elements that are proportional to pure, single-mode 
Gaussian states (i.e., rank-one Gaussian operators) with covariance matrix 


P-B 


cos 9 b 
sin 9 b 


— sin 9 b 
cos 9 b 


Lb 0 A 

0 ^/LbJ 


cos 9 b 
— sin 9 b 


sin6*s \ 

cos 9 b I ’ 


(5) 


which is, in fact, the covariance matrix of a squeezed-vacuum state. The various outcomes 
of the measurement correspond to the points {6} = {xb,Pb} in the phase plane; the 
corresponding POVM elements are obtained by displacing the squeezed-vacuum state to all 
points in the phase plane. (More generally, a single-mode Gaussian measurement can have 
POVM elements that are Gaussian convex combinations of the rank-one POVM elements, 
i.e., that are proportional to mixed, single-mode Gaussian states, but such measurements are 
noisier versions of the rank-one Gaussian measurements and thus are never optimal for our 
considerations.) Homodyne measurement has Lb = O', heterodyne measurement has = 1; 
and for measurements in between, 0 < Lb < 1. 

After performing such a measurement on B with outcomes b, the conditional state pA\b 
has mean quadratures that depend on b, but its covariance matrix, given by 


(TA = A-CiB + p,B)-^C^, 


( 6 ) 


is independent of b ll25ll . Thus, the eigenstates of pA\b are, in general, displaced squeezed 
number states; measuring in this basis minimizes the the Shannon entropy of the outcome 
probability distribution, making it equal to the von Neumann entropy S{A\b) of the 
conditional state pA\b- This von Neumann entropy is given by S{A\b) = F(-\/det cta) BUI, 
with 


F{x) 


X -1-1 
2 


X + 1 
2 



(7) 


The Gaussian quantum discord (GQD) is then given by 


Dgqd(5 ^ A) = H^jL^^\A\B) - 5(A|5), (8) 


where the classical conditional entropy, 

^|min)(^|^) ^ niinS'(A|6) = min f(s/ det cta ), 

Ub Ub \ / 


( 9 ) 
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is now obtained by minimizing S'( 74 | 6 ) only over Gaussian measurements on B. 

For Gaussian states in the standard form Q and having a = 6 = c+ l, itis interesting 
to note that the quantum conditional entropy of the state with d = c, referred to as the 
correlated-correlated (CC) state, is smaller than the quantum conditional entropy of the state 
with d = —c, referred to as the correlated-anticorrelated (CA) state. On the other hand, the 
classical conditional entropy @ is the same for these two separable states ffTTl . This implies 
that the GQD of the CC state is larger than that of the CA state, although the marginal states 
of these two separable states are the same. Note that given c, the CC and CA states are the 
nonentangled states that have maximal correlations in their quadratures. 

Recently, using a connection between the continuous (differential) Shannon entropy 
of the Wigner function and the Renyi-2 entropy, Gaussian Renyi-2 discord was defined as 
a measure of nonclassical correlations in Gaussian states (211 . In this measure, the von 
Neumann entropies in Eq. ([^ are replaced by Renyi-2 entropies 

D 2 {B ^ A) = min S 2 {A\b)-S 2 {A\B), (10) 

fJ'B 

where S' 2 (A| 6 ) = —ln(Tr[p^|J) = ^ ln(det cTyi) and S' 2 (A|i?) = | ln(det cr^^/det B). 
Notice that, the conditional entropy S 2 {A\b) corresponds to the continuous Shannon entropy 
of the Wigner function of pA\b up to a constant (24l . There is, however, no Gaussian 
measurement whose outcome probability distribution is equal to the Wigner function, as 
the noncommuting observables xa and pA cannot be measured simultaneously without some 
noise penalty. The CC and CA states have the same Gaussian Renyi-2 discord, and Gaussian 
states with no correlations in one of the quadratures (d = 0) have zero Gaussian Renyi-2 
discord. 

Neither the GQD nor the Gaussian Renyi-2 discord satisfy the condition of nonclassical 
correlations lfT3l for Gaussian protocols, because they use a nonGaussian measurement on A. 
We turn now to formulating an operational, discord-type measure of nonclassical correlations 
for Gaussian states that is based purely on Gaussian measurements. 

3. Operational Gaussian discord 

We refer to our new measure as operational Gaussian discord (OGD) and define it as 

Dogd{B ^A)= Hifr\A\B) - Hi^^\A\B), (11) 

where {^\B) is the minimum conditional entropy of A after performing local Gaussian 
measurements on A and B, and H^j'^\A\B) is the minimum conditional entropy of the same 
subsystem after performing a joint Gaussian measurement on A and B. The entropies are 
the continuous (differential) Shannon entropy of Gaussian probability distributions, which 
for a single mode are given by | ln(det a) + ln(27re), with a being the covariance matrix 
of the probability distribution (see Appendix A). In our notation, A and B denote that the 
entropies are calculated using outcome probability distributions of the measurements. As all 
the probability distributions are Gaussian, in order to calculate the OGD ([TT]), one just needs to 
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minimize the determinants of the eovariance matriees of the eonditional Gaussian probability 
distributions for the outcomes of local and joint Gaussian measurements. For a discussion 
of Gaussian measurements, conditional probability distributions, and the corresponding 
entropies, see Appendix A. 

In general, the POVM elements of a two-mode Gaussian measurement are proportional 
to two-mode Gaussian states whose covariance matrix, according to the Williamson 
theorem [l27ll . can be written as 

Pj = 021)8. (12) 


Here the 2x2 identity matrix 1 represents the single-mode vacuum state (the choice of units 
can be thought of as setting h = 2), i/il © Z/ 2 I corresponds to product thermal states with 
variances ui and 02 in the two modes, and S represents a symplectic transformation [l28ll . For 
our purpose, that is to minimize the entropies of the outcome probability distributions, we 
consider Gaussian measurements that have rank-one POVM elements, i.e., oi = 02 = 1; 
thus, the POVM elements are proportional to pure Gaussian states. Measurements with 
mixed POVMs will add more noise and increase the entropy. Any symplectic matrix can be 
expressed as S = K[s(ri) © s(r 2 )]L, where K and L represent beamsplitter transformations 
and s(r) represents a single-mode squeezing operation [[291 Using this expression for 


S in Eq. (12) and knowing that the action of a beamsplitter on vacuum states results in 
vacuum states, we can write the covariance matrix of the POVM elements of a two-mode, 
joint Gaussian measurement in the form 


0b)B'^( 77)(^^ © AiB)B(?7)R(0A, fs) 




(13) 


where 


B{ri) 


0 

\/l 


0 

Vv 

0 


\ 0 

describes a beamsplitter transformation. 


-v^T^ 0 \ 

0 -s/T^ 

Vv 0 

0 Vv / 


4>b) 


^COS (pA 
sin cpA 
0 

V 0 


— sin (pA 
cos CpA 


0 

0 


0 

0 \ 

0 

0 

cos (pB 

— sin (pB 

sin (pB 

cos (pB / 


(14) 


(15) 


describes pre-beamsplitter single-mode phase shifts, and is defined as in Eq. ([^, with p,A 
defined analogously. As we see from the above expression, the joint Gaussian measurement 
corresponding to this covariance matrix can be realized by two phase shifters and a 
beamsplitter followed by two local Gaussian measurements. Obviously, for cpA = (pB = 0 
and ?7 = 1, we obtain the covariance matrix of a local Gaussian measurement. 


A^l = A^a © A^s- 


( 16 ) 
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As shown in the appendix, after performing a joint Gaussian measurement, the 
eovarianee matrix of the eonditional probability distribution for A is given by 

dTA^j = A - (17) 


which is obtained from a joint Gaussian probability distribution with the covariance matrix 


AB,J — (J'AB + P'J 


I A + C + j 

+ Hqj B + pb,j j 



(18) 


After local Gaussian measurements on A and B, the covariance matrix of the conditional 
probability distribution for A is 


(Ta,l = A + Pa- C(B + pb) ^C'^. 


(19) 


Thus, by using Eqs. ( [T9l ) and ( [T7| ), the OGD measure becomes 

Dogd{B -a- A) = min ^ In (det ct^.l) - min ]- In (det cr^,j). (20) 

ua,ub Z uj Z 

Operational Gaussian discord is always nonnegative, because the set of all joint 
measurements includes all local measurements; hence, the conditional entropy minimized 
over all possible joint measurements, H^f^AlB), can never be larger than the conditional 
entropy miminized over all local measurements, i7“]f(A|S). In addition, the OGD of 
product states is zero. In this case, for a joint measurement we have = A + pf with 
p'a = Pa,j — + pb,j)~^Pcj, which is equivalent to a local measurement on A with 

covariance matrix p'^', hence, det dr^ j cannot be smaller than det a-A,L- 

4. Examples: operational Gaussian discord for some Gaussian states 

In general, it is not clear how to obtain a closed-form expression for OGD for 
an arbitrary Gaussian quantum state. The conditional entropy with local Gaussian 
measurements must be minimized over four parameters, {6a, Ob, La, Lb}, and the 
conditional entropy with joint Gaussian measurements must be minimized over seven 
parameters, {(pA, 4>B,f], 6 a, 6b, La, Lb}- In the following, by using analytical and numerical 
methods, we calculate OGD for some Gaussian states in the standard form Q. 

4.1. Entangled and separable Gaussian states 

Let us consider a class of Gaussian states whose covariance matrices Q are parameterized by 
a and t such that a = b > 1 and c = —d = ts/a? — 1 > 0, where 0 < f < 1. When t = 1, this 
is a pure two-mode squeezed-vacuum state, with a = b = cosh2r and c = —d = sinh2r, r 
being the squeezing parameter. For t > ■\/{a — l)/{a + 1), the state is entangled ll23ll : for the 
other values of t, the state is separable. The boundary between separability and entanglement, 
i.e., t = \/{a- l)/(a -f 1), is occupied by the CA state. 
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By minimizing det a-A,L over all local Gaussian measurements we find that the optimal 
local measurements for all values of a and t are two heterodyne measurements, i.e., La = 
Lb = I- This gives a symmetric covariance matrix for the conditional probability distribution: 


(^A,L 



a)f^ 0 

1 + a + (1 



( 21 ) 


In order to minimize detdr^ j one can guess that, as the quadratures are equally 
correlated but with a different sign, the covariance matrix of the POVM elements of the joint 
Gaussian measurement must be in the same form, with the sum ctab + I^j = ^ab,j enhancing 
the correlations and minimizing the determinant of the conditional covariance matrix j. 
This means that p, j must be the covariance matrix of a two-mode squeezed state, as it is the 
only pure Gaussian state in that form. 


fJ'AJ — l^B,J 


- 2 


Ul/L + L) 


Ul/L + L) 


l^c,j = 


- 2 


\{l/L-L) 


-\{l/L-L) 


( 22 ) 

(23) 


i.e., (t)A = (t>B = Oa = 0, 6b = 7r/2, p = 0.5, and La = Lb = L. Numerical calculations 
confirm that this is the optimal choice for the joint Gaussian measurement. Minimizing 
det d-A,j over the parameter L gives 


L 


{ 1 — — a — -f 2ts/a? — 1 

1 + t‘^ — a {1 — t'^) 

0 . 


0 < t < s/{a - l)/(a -f 1), 
^/{a — l)/(a -f 1) < f < 1. 


(24) 


The corresponding covariance matrix of the conditional probability distribution for 0 < t < 

A/(a - l)/(a -f 1) is 


(1 -f a) (1 — f^) 0 

* 0 (1 + a) (1 - f^) 


and for other values of t is 

^A,J = 2 


a — ty/a^ — 1 0 

0 a — tyj — I 


(25) 


(26) 


For the CA state and entangled states, an optimal measurement is a beamsplitter followed 
by two homodyne measurements (POVM elements are two-mode infinitely squeezed states). 
For the separable states, the local measurements after the beamsplitter are measurements in a 
displaced squeezed-vacuum basis, varying between heterodyne for t = 0 and homodyne for 
t = y/ {u — l)/(fl + !)• 
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Figure 1. The solid blue line shows the operational Gaussian discord (OGD), and the dashed 
red line shows Gaussian quantum discord (GQD) for a Gaussian state in the standard form Q, 
with a = b = 10 and c = —d = ty/OO. The optimal local Gaussian measurements are two 
heterodyne measurements for all values of t. For t > -^9/11, the optimal joint Gaussian 
measurement is a 50:50 beamsplitter followed by two homodyne measurements. Notice that 
for t > -^9/11 (vertical line), the state is entangled. For other values of t, the optimal 
joint Gaussian measurement is local measurements in a displaced squeezed-vacuum basis after 
the beamsplitter, varying between no squeezing (heterodyne) for f = 0 to infinite squeezing 
(homodyne) for t = ^JO/11. 


The OGD for these states is given by 


DoGoi^B —>■ ^4) — ln(l T Cl + — cit^) 

_ |ln((l + a)(l 

1 ln(2a — 2t\/a^ — 1), 


0 < t < \/{a- l)/(a + 1), 
<t <1. 


(27) 


For these states Dgqd{B A) < Dogd{B —)■ ^4), which we illustrate in a particular 
case in Fig. This implies that the difference between the classical and quantum conditional 
entropies is less than or equal to the difference between the conditional entropies obtained by 
local and joint Gaussian measurements. Also, for the two-mode squeezed-vacuum state, we 
have DoGoiB A) = 2r; the quantum discord of this state is equal to the von Neumann 
entropy of the marginal state, S{B), which for large values of r scales as 2r + 1 — 2 In 2. 


4.2. Correlated-correlated and correlated-anticorrelated states 

Here we consider separable Gaussian states parameterized by c and q such that a = 6 = c + 1 
and d = qc, with c > 0 and — 1 < g < 1. The parameter q controls the correlation in the 
p-quadratures; by changing q from —1 to 1, the state changes from the CA state to the CC 
state. 

We first minimize det a-Ap for the minimum conditional entropy with a local Gaussian 
measurement. Numerical calculations show that 6a = Ob = 0, as expected because the state 
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is in standard form. The minimizing values of La and Lb are 


La = 


2 + c — cq^ 


and 


Lb — 


(l + c)(g^-l) + c|g|V4 + 2c-2cg^ fl + < lol < 1 

(l + c)2-(l + c2)g2 ’ + <191 Si, 

0 < |g| < (l + 2 c)-i/S 


(28) 


(29) 


0 , 


and these give 


det (Ta,l = \ l + La + 


c(l + Lb) 
1 + c + T 


B 


1 + C+ — 
La 


c^qM 


B 


(1 + c)Lb + 1 


(30) 


Aeeording to the above expressions, the optimal loeal measurement on A for all values of 
g is in a displaeed squeezed-vaeuum basis, whieh limits to a heterodyne measurement when 
|g| = 1. The optimal loeal measurement on i? for (1 + 2c)“^/^ < |g| < 1 is also in a displaeed 
squeezed-vacuum basis, but for the small eorrelations in the p-quadratures, |g| < (l + 2c)“^/^, 
the optimal local measurement is homodyne. For |g| = 1 the local measurements on both A 
and B are heterodyne measurements; in this case, we have = 4(1 + c)/(2 + c). 

Using numerical calculations, we find that the POVM elements of the optimal joint 
Gaussian measurement are two-mode squeezed states, with covariance matrix p,j given by 


Eqs. ( |22| ) and ( [23] ). We obtain 

4(1 L)[l L 2cL)[l L cL — cqL)[l c L cq) 
=- (1 + 2L + 2cL + ^ 


which is minimized by 

L = ■?-1+V4 + 4c-4c^^ 

3 -f 2c(l -q) - q 

The expression for L shows that for g = 1 we have L = 1; i.e., the optimal joint Gaussian 
measurement is a 50:50 beamsplitter followed by two heterodyne measurements. In this case, 
Pj = 1 <3) 1, and this measurement is equivalent to two local heterodyne measurements. 
Moreover, it is easy to see that for (1 -f 2c) < g < 1, the minimum conditional entropies 

with local and joint Gaussian measurements are the same, det o-aj = det and thus OGD 
is zero. We also observe that the parameter L decreases as g decreases. For g = —1 we have 
L = 0 and y/detd^^ = 2, which corresponds to performing two homodyne measurements, 
with 9a = 0 and 9b = j‘I, after the beamsplitter. 

In Fig. [^ we compare GQD and OGD measures for Gaussian states with c = 9. As 
shown, GQD for the state with g = 1 is larger than for the state with g = —1. In contrast, 
according to OGD, the state with g = 1 has zero correlation, and the state with g = — 1 has 
the maximum correlation. 
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Figure 2. Operational Gaussian discord (OGD) (solid blue line) and Gaussian quantum 
discord (GQD) (dashed red line) for Gaussian states with a = 5 = c+ l = 10 and d = Qg. 
The parameter q controls the correlation between the p-quadratures of the joint system. OGD 
monotonically decreases in the interval —l<q< l/v/l9, and for l/v/l9 < g < 1 OGD is 
zero; note that q = l/v/l9 is the point at which the optimal local measurement on B changes 
to homodyne measurement. According to GQD, the CC state, q = 1, has more nonclassical 
correlation than the CA state, q = —1, but the OGD measure attributes zero nonclassical 
correlation to the CC state and the maximal nonclassical correlation in this class to the CA 
state. 


4.3. Asymmetric Gaussian states 

Bipartite Gaussian states whose marginal states are not the same are asymmetric. To explore 
properties of such states, we consider separable, asymmetric Gaussian states in the standard 
form, which are parametrized by b, v, and s with a = b + v, c = |s|, and d = s, where b > 1, 
n > 0 and c = |s| < 6 — 1. 

Using numerical and analytical calculations, we find that the optimal local measurements 
that minimize det a-A,L for all values of s and v are heterodyne measurements. La = Lb = 1, 
which yields a symmetric covariance matrix for the conditional probability distribution. 



(33) 


In order to minimize det cta,./, we consider the cases s > 0 (c = d) and s < 0 (c = —d) 
separately. Note that for s = 0 the state is a product state. For c = d, as for the CC 
state in the previous subsection, we find that for all values of |s| and v the optimal joint 
Gaussian measurement is a 50:50 beamsplitter followed by two heterodyne measurements 
if A = (j)B = Oa = Ob = 0, La = Lb = I, V = 1/2), which implies j = For the 
case c = —d, however, we find that the optimal joint Gaussian measurement is described by 
parameters ^a = fs = Oa = 0, 9b = vr/2, rj = 1/2, and La = Lb = (b — 1 + s)/(b — 1 — s); 
i.e., the POVM elements are two-mode squeezed states. In this case, we obtain a symmetric 
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covariance matrix for the eonditional probability distribution as 

_ + b + V — s^/(b — 1) 0 

y 0 1 + b + V — s^/(b — 1) 

As a eonsequence, OGD for these states is given by 


(34) 


-Dogd(-B A) 



In 1 + 


(1 + b)(b'^ — — 1 — V + bv) J ’ 


.2 


0 < s < 6 — 1. 


1 — 6 < s < 0. 


(35) 


Notiee that the optimal loeal and joint Gaussian measurement strategies are independent of 
the value of v, and for n —)■ cx), OGD is zero. 

5. Operational significance 

We now present the operational significance of our measure in terms of a Gaussian protocol 
for encoding information onto Gaussian quantum states (see Fig. 1). In this protocol, two 
independent classieal random variables, Xg and Ps, represented by the veetor Xg = (xg,pg) 
and deseribed by Gaussian probability distributions with the same varianee Vg, are eneoded 
on the X- and p-quadratures of subsystem A of a joint system in the Gaussian state pab- The 
eneoding proeedure is done by applying the displaeement operator DA{Xg) = eyi'p[i{pgXA — 
XsPa)/2] and averaging over the Gaussian distributions for Xg and pg. The state after eneoding 
thus becomes 



(36) 


The state pf^ is also Gaussian, with eovarianee matrix = ctab + Vgl © 0, where 0 is 
the 2x2 zero matrix. 

The aim is to obtain an estimate of the signals, = {xe,Pe), by using some 
measurement strategy that takes advantage of the correlations between the subsystems in 
such a way that the elassical mutual information I{Xg^Yf) is maximized. It was shown, 
using Holevo’s theorem, that for maximal eneoding, i.e., K —)■ cxd, the differenee between the 
maximum extraetable information with and without restrieting to loeal measurements is equal 
to quantum diseord of the state pab dill- In order to saturate the extraetable information, 
however, nonGaussian measurements are required, in the way we described earlier for 
quantum diseord. Here we eonsider a Gaussian version of the protoeol, in which there are 
two measurement strategies: loeal Gaussian measurements and joint Gaussian measurements. 

Consider first the ease where subsystem B is not available to us. Assuming the state pab 
was in the standard form, the marginal state pA is symmetric with variances a in the x- and 
p-quadratures. In this case, the eovarianee matrix of the outcome probability distribution of a 
Gaussian measurement is given by 



( 37 ) 
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Figure 3. By applying a Gaussianly distributed local displacement operator, a pair of classical 
Gaussian random variables Xg = (xs,Ps), both having the same variance Vg, are encoded on 
X- and p-quadratures of subsystem A, which is part of a bipartite system in Gaussian state with 
the covariance matrix (tab El- (a) In the first strategy, one performs optimal local Gaussian 
measurements, whose POVM elements are described by the covariance matrix pA 0 Ms- 
After post-processing the data, one obtains a signal estimate Yf, = {xe,Pe) such that the 
mutual information Il{Xs, Yf) is maximized. The covariance matrix of the joint probability 
distribution is = (tab + Pa (B Pb- (b) In the second strategy, one performs an optimal 

joint Gaussian measurement such that the mutual information Ij{Xs, Yf is maximized. As 
shown in the text, the most general form of a joint Gaussian measurement consists of two 
phase shifters and a beamsplitter (BS) followed by two local Gaussian measurements. The 
covariance matrix of the probability distribution for the outcomes of this measurement is 
^'ab j ~ ^AB + MJ’ where pj is the covariance matrix of the POVM elements. In the 
limit of maximal encoding (14 -A oo), the difference between Ij{Xg, Yf and G{Xg, Yf is 
equal to operational Gaussian discord (OGD) of the state pab- 


By using the expression for the mutual information of two parallel Gaussian ehannels [l3T]| . 
the mutual information between Xg and an estimate of it, Y^, that is obtained after the 
measurement is given by 


I{Xg,Yf 


-In 1 + 


14 


a, + La 


+ 




14 


0 + 1 / La 


(38) 


This quantity is maximized for La = 1, i.e., by performing heterodyne measurement. There 
are two sourees of noise redueing the mutual information: the noise of the quantum state and 
the noise assoeiated with the measurement. While the former noise is inevitable due to the 
uneertainty prineiple, the latter eould be redueed if subsystem B was available to us. In that 
ease, one eould take advantage of the eorrelations by performing some measurement on B 
and post-proeessing the outeomes in order to effeetively reduee the noise of the state, thus 
allowing extraetion of more information about the signals. 

When both subsystems are available, in the first strategy, one performs loeal Gaussian 
measurements on subsystems B and A. This yields a eonditional probability distribution for 
A, with eovarianee matrix 0-^,1 = ^a,l + 141, where d-A,L is given by Eq. (191. The mutual 
information, given by/(X^, 14 ) = l[ln(l + 14 /api)+ ln(l+ 14 /a£, 2 )], where opi andap 2 are 
the eigenvalues of should be maximized over the loeal measurements, whieh means to 
maximize it over the loeal-measurement eovarianee matriees p,A and ps- Thus, the quantity 
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of interest is 


hiXs^Ye) = max - 


In ( 1 + —^ + In (^1 + — 

ail) V Q>L2 


(39) 


In the second strategy, by using a joint Gaussian measurement, one obtains a conditional 
probability distribution with covariance matrix af j = j + Vgl, where o-aj is given by 
Eq. (17). The mutual information, I{Xs,Ye) = ^[ln(l + Vg/aji) + ln(l + Vs/aj 2 )], where 
aji and aj 2 are the eigenvalues of cr^ j, should be maximized over the covariance matrix p,j 
that describes the joint Gaussian measurement, so the quantity of interest is 


Ij{Xs,Ye) = max- 

M.7 z 


In ( 1 + —^ + In (^1 + — 

Ojl ) V 


(40) 


In the limit that the classical signals have very large power, Vg is much larger than any of 
the eigenvalues. In this situation, we have 


lL{Xg,Ye) ~ InK - min ^^(0^10^2) = InK - min ^ ln(det cr^,z,) 
Ua,Ub Z Ua,Ub Z 


(41) 


for the first strategy, and 


Ij{Xg, Ye) ~ In 14 - min ^ ln(ajiaj2) = In 14 - min ^ ln(det (Ta,j) 

uj Z uj Z 


(42) 


for the second strategy. In the limit I 4 —)■ C) 0 , the difference between these two mutual 
informations is equal to the OGD of pab. 


Ij{Xg, Ye) - hiXg, Ye) = DoGuiB ^ 4). 


(43) 


This relation provides the operational significance for our measure. 

For some Gaussian states the local and joint Gaussian measurements used to minimize 
the conditional entropies for the OGD of are independent of the value 14, as shown for 
the states considered in Sec. |4.3[ whose conditional probability distributions are symmetric 
Gaussian functions. In this case, one can easily see, for example, by considering Eqs. 
and (|34]), that 


Ij{Xg, Ye) - hiXg, Ye) = DocDiB ^ 4) - ^ (44) 

where —)■ A), the OGD for the state encoding, is zero for maximal 

encoding. According to this relation, the difference between mutual informations obtained 
by the above strategies is equal to the amount of nonclassical correlation in terms of OGD 
consumed by encoding the signal. This implies that, for any value of 14, there is no difference 
between the two strategies for the CC state, IiiXg, Ye) = Ij{Xs, Ye); however, for the CA 
state the joint Gaussian measurements strategy is always advantageous with respect to the 
local Gaussian measurements strategy, Il{Xs, Ye) < Ij{Xs, Ye). 
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6. Conclusion 

We have shown that operational Gaussian diseord (OGD) is a new diseord-type measure of 
nonelassieal eorrelations for Gaussian states, whieh ean be experimentally measured by using 
local and joint Gaussian measurements. We have demonstrated an operational significance 
for this measure in terms of a Gaussian quantum protocol for extracting information about 
a classical signal encoded on one subsystem; for maximal encoding, OGD is the additional 
accessible information that comes available when an experimentalist throws off the shackles 
of local Gaussian measurements and starts using joint Gaussian measurements. This measure 
might also be useful for quantifying nonelassieal correlations in resources of other Gaussian 
protocols that involve Gaussian states, Gaussian operations, and Gaussian measurements. 

An interesting open question is how to define a similar measure for discrete-variable 
systems. This measure can be defined as the difference between two conditional entropies of 
one subsystem minimized by local and joint measurements. Such a measure might have an 
operational significance in terms of the discrete-variable version of the protocol we considered 
in this paper and other quantum protocols. 
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Appendix. Gaussian measurements and entropy of Gaussian probability distributions 

In quantum mechanics, the uncertainty principle only allows noisy simultaneous 
measurements of the phase-space quadratures [|^ . In the general formalism for phase-space 
measurements, the POVM elements associated with the outcomes (x,p) of a single-mode 
measurement are given by 

n(a:,p) = ^D{x,p)TIqD\x,p), (.1) 

where D{x,p) is the displacement operator, J dx dpll{x,p) = 1, and the quantum state Ho, 
which can be assumed to have zero first-order moments, is a characteristic of the measurement 
device, sometimes called a quantum filter or ruler [|33l|34l. The phase-space measurements 
for which Ho is a Gaussian state are called Gaussian measurements Il28l : we show in the 
following that the outcome probability distributions are Gaussian. It was shown, in general, 
that Gaussian operations can be implemented using linear-optical elements and homodyne 
measurements [|^ |2^ . Hence, Gaussian measurements are equivalently defined as a set 
of measurements that can be implemented using Gaussian ancilla states, Gaussian unitary 
operations, and homodyne measurements [l35l[36]l . As discussed in Sec.|^ we are interested 
in rank-one Gaussian measurements, i.e.. Ho is a squeezed-vacuum state, as these states satisfy 
the minimum uncertainty relation, and the measurement is as accurate as possible. 
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Note that Gaussian measurements can be put in a slightly more general framework that 
uses measurements with noncovariant POVMs Y\.{x,p,y) = p{y)D{x,p)Y\.Q{y)D''{x,p)/ATT, 
where the Gaussian states no(r/) are labeled by an additional parameter y, and p{y) is the 
probability of outcome y, since / dx dp n(x, p, y) = p{y)l ifTSll . These measurements can be 
thought of as performing one of a set of random Gaussian measurements by flipping a coin, 
governed by the probability p{y), to choose which Gaussian measurement. Since we optimize 
over all Gaussian measurements, a random choice of Gaussian measurements would not be 
optimal. Formally, including the parameter y in the minimization of the entropy would pick 
out the value of y for the optimal measurement and lead to the same minimum entropy. 

As discussed in the main text, a two-mode Gaussian measurement can be implemented 
using linear-optical elements and single-mode Gaussian measurements; the POVM elements 
are 

nCJC) = (.2) 


Here Ho is a Gaussian state with zero mean quadratures and the covariance matrix pLj of 
Eq. (13); X = {X) = {xa,Pa,xb,Pb) is the vector of mean quadratures of the state n(X) 
and represents the outcomes of the measurement; ^ = (^i, ^ 2 , ^ 3 , ^ 4 ) G and ^ = —XJ/2 
with 


J = 




/ 0 

1 

0 

o\ 

[ Ja 

0 A 

-1 

0 

0 

0 


1 

0 

0 

0 

1 



1 0 

0 

-1 

0/ 


(.3) 


being the fundamental symplectic matrix; and 


D{$,) = 


- = ^-iXJXTl2 _ 


= e 


i{PAS^A-^APA+PBS^B-XBPB)/‘i _ 


= D{X) 


(.4) 


is the two-mode displacement operator. 

By using POVM elements ([^ The probability distribution of the outcomes of the 
Gaussian measurement performed on a bipartite quantum system in state pab is given by 


P{X)=Tt[pabII{X)]. 


(.5) 


By using the characteristic function of the state, ^ab{^) = T''<^[pabD{^)], and the 

characteristic function of the POVM-element states, <hn(x)(^) = exp(—-f 

i^X^) /(IGvr^), this distribution can be written as 

nx) = (.6) 


For a zero-mean Gaussian state pab with covariance matrix ctab, the characteristic function 
becomes = exp(—so we have 


P{X) 


1 




dvr^-^/det a-AB 


(-7) 
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Thus the probability distribution is a Gaussian function with covariance matrix (Tab,j = 
f^AB + as in Eq. ( [T^ . 

Using the continuous Shannon (differential) entropy, the entropy of the joint probability 
distribution can be found as 


H{A,B) = - / d^X P{X)\n{P{X)) = - ln(det + 21n(27re). 


(. 8 ) 


The constant 2 ln(27re) does not have an absolute significance; the continuous entropy is only 
defined up to an additive constant. The difference between two such entropies does, however, 
have an absolute significance. 

When the covariance matrix o-abj is written in terms of the block matrices A, B, and 
C of Eq. (181, the inverse is given by lUTlI 




—cr 


^A,J 

C^A-i 


where 


B,J' 


<^A,J 


^b]j 


^A,J 



= A-CB 


= B - C^A-^C 


(.9) 


(. 10 ) 


By using this expression and the probability distribution 0. one can easily find the 
conditional probability distribution for A, given measurement results on 5, as 


P{xa,Pa\xb,Pb) = P{Xa\Xb) = 


e 2 " 


(. 11 ) 


27ry^det cta . j ' 

where R = Xa — XbB^^C^. The covariance matrix of this distribution, cta.j, is 
independent of the outcomes xb and pb- Hence, the continuous Shannon entropy of the 
conditional probability distribution is given by 


H{A\B) = - ln(det (Ta,j) + ln(27re). 


The conditional entropy (.12) can also be written as 


H(A\B) = fn 


det crAB,j 
det B 


+ ln(27re), 


(. 12 ) 


(.13) 


since for joint classical probability distributions we have H{A\B) = H(A, B) — H{B). This 
relation can be regarded as a consequence of the identities 


det (Tab,, 
det B 


= det A I 1 - Tr[CB-^C^A-^] + 


(det C)^ 
det B 


= det (Ta,j. 


(.14) 


We note that by setting pj = 0, i.e., by removing the tildes on all quantities so that 
^AB,j = <^AB, the probability distribution P{X) of Eq. ([^ becomes the Wigner function 
of the joint state pab, with Eq. ([^ giving the continuous Shannon entropy of the Wigner 
function. Moreover, we can model a local measurement on B by removing the tildes from 
A and C, i.e., by setting pa,j = 0 = pc,j so that cta,j = (Ta [see Eq. ([^]; in this case the 
conditional distribution ( |.11[ ) becomes the conditional Wigner function of A, with Eq. (.12) 
giving the corresponding continuous Shannon entropy. 
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